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Abstract We present a study of the relationship between galaxy colour, stellar mass, 
and local galaxy density in a deep near-infrared imaging survey up to a redshift 
of z ~ 3 using the GOODS NICMOS Survey (GNS). The GNS is a very deep, 
near-infrared Hubble Space Telescope survey imaging a total of 45 arcmin 2 in the 
GOODS fields, reaching a stellar mass completeness limit of = 10 9 ' 5 M© at z = 
3. Using this data we measure galaxy local densities based on galaxy counts within 
a fixed aperture, as well as the distance to the 3"', 5 and 7 th nearest neighbour. We 
find a strong correlation between colour and stellar mass at all redshifts up to z <~ 3. 
We do not find a strong correlation between colour and local density, however, the 
highest overdensities might be populated by a higher fraction of blue galaxies than 
average or underdense areas, indicating a possible reversal of the colour-density 
relation at high redshift. Our data suggests that the possible higher blue fraction at 
extreme overdensities might be due to a lack of massive red galaxies at the highest 
local densities. 



1 Introduction 

Despite a wealth of recent studies, disentangling the influence of stellar mass and 
environment on galaxy formation and evolution remains a highly debated topic. The 
presence of a strong correlation between galaxy colour and stellar mass is supported 
by numerous observational studies in the local universe (see e.g. lfl"2l fPJI ) and up 
to intermediate redshifts of z ~ 1 (see e.g. (5] [To)). However, galaxy colours not 
only depend on mass, but also on a galaxy's environment. The preference of red 
galaxies for denser local environments, first noticed by fl31 . and confirmed by (8), 
is now well-studied in the local universe (e.g. |[T4l[T8l ). In the early universe, how- 
ever, the evidence is controversial. While some studies find that the colour-density 
relation at z ~ 1 is mainly due to a bias in stellar mass selection and only persists 
for low-mass galaxies (e.g. lfT7l[TTl[T0l ). others argue that a strong colour-density 
relation is already in place at z ~ 1.4 (e.g. |5] [9)) even at fixed stellar mass (5|. 
In this study we investigating the relationship between galaxy colours, stellar mass 
and local densities in the critical redshift range at z > 1.5, probing the period in 
which most of galaxy formation takes place. We use data from the GOODS NIC- 
MOS Survey (GNS), a large Hubble Space Telescope survey reaching a stellar mass 
completeness limit of ~ 10 95 M at z ~ 3. 
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2 Data and Analysis 
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In this study we use data from the GOODS NICMOS Survey (GNS) covering a 
total area of about 45 arcmin 2 in the GOODS fields (4}. The limiting magnitude 
reached at 5(7 is Has = 26.8, which is more than 2 mag deeper than ground based 
near-infrared imaging lfl6l . To obtain photometric redshifts, stellar masses and rest- 
frame colours, the NICMOS H-band sources were matched to a catalogue of optical 
sources in the GOODS ACS fields, yielding a BVizH photometric catalogue. 

The photometric redshifts are measured with the HYPERZ code (T). F° r the 
high redshift sample (z > 1.5) we obtain a photometric redshift uncertainty of 
Gaz/(\+z) = 0-10' while galaxies atz < 1.5 show a slightly lower, but still compara- 
ble scatter of Oaz/(i+z) = 0.08. The photometric redshift uncertainties are used in a 
set of 100 Monte Carlo simulations to estimate the uncertainties in local density. For 
this purpose we randomize the photometric redshift input according to the photo-z 
error, accounting for scattering in and out of the redshift range, as well as dealing 
separately with catastrophic outliers. The randomized photo-zs are then used to re- 
compute the local densities, giving typical average uncertainties of ~ 0.2 dex. We 
plot the results of the simulations in Figures[2]and[3]discussed in Section [3~2l 

The stellar masses and (U — B) rest-frame colours are measured by fitting a large 
set of synthetic spectra constructed from the stellar population models of (2), as- 
suming a Salpeter initial mass function. The average uncertainty of our stellar mass 
estimate is ~ 0.3 dex. 

We compute local densities using the fixed aperture method as well as the N th 
nearest neighbour approach. The absolute densities are corrected for edge effects 
and normalized by the median value in a redshift bin of Az = 0.25, yielding a relative 
local density log (1 + 5n)- Comparing the different density estimators we conclude 
that the 3 nearest neighbour density has the largest dynamical range and only 
slightly higher uncertainties and is therefore best suited to trace the extremes in the 
local density distribution. We show the results for log (1 + 83) in the following. 



3 Results 

3.1 Colour-magnitude relation 

Figure Q] shows the (U —B) - Mb colour-magnitude relation in three different red- 
shift bins: 1.5 < z < 2, 2 < z < 2.5 and 2.5 < z < 3. In the three top panels the 
sample is split in 4 stellar mass bins with a bin size of 0.5 dex; the symbols are 
shaped and sized according to the galaxy's stellar mass: larger symbols represent 
higher stellar mass galaxies. In the three lower panels the sample is split in bins of 
local density with a bin size of 0.6 dex. Here, larger symbols correspond to galaxies 
located in higher overdensities. The red solid lines indicate the expected location of 
the red sequence observed at z ~ 1, and evolved passively back in time according 
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Fig. 1 Colour-magnitude relation in three redshift ranges, divided in bins of stellar mass (top 
panels) and local density (bottom panels). The larger symbols correspond to higher stellar mass 
in the top panels, and to higher overdensity in the bottom panels. The expected location of the 
red sequence at the respective redshift is shown as red solid line, while the limit to the blue cloud 
below which a galaxy is considered blue is plotted as blue long-dashed line. 



to the models of |fl9l . The blue long-dashed lines is the border to the blue cloud, 
defined as being 0.25 magnitudes bluer than the red sequence. The different slices 
in stellar mass are clearly separated in colour-magnitude space. The separation is 
perpendicular to the red sequence, and the distance from the expected red sequence 
increases with decreasing stellar mass. The different slices in local density on the 
other hand are not well separated in colour-magnitude space. The relation between 
colour and local density is much weaker than the relation between colour and stellar 
mass. 



3.2 Colour-density relation and the role of stellar mass 

The left panel of Figure|2l shows rest-frame (U —B) colour (top) and the fractions 
of blue and red galaxies (bottom) as a function of stellar mass in three redshift bins. 
The correlation between (U —B) colour and stellar mass log M t is highly significant 
at all redshifts. The cross-over mass, i.e. the stellar mass at which there is an equal 
number of red and blue galaxies, stays roughly constant at log M* ~ 10.8. This is 
the same cross-over mass found at lower redshifts of 0.4 < z < 1 fTOl . 

The right panel of Figure [2] shows the relation between (U — B) colour (top row) 
and the fraction of blue and red galaxies (bottom row) as a function of relative over- 
density log (1 + In this figure we plot the results of the Monte Carlo simulations 
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Fig. 2 (U — B) colour and fraction of red and blue galaxies as a function of stellar mass (left) 
and local density log (1 + 53) (right). Each line corresponds to the results of one Monte Carlo 
Simulation, averaged in bins of stellar mass (left) and local density (right). 



rather than the original datapoints. Each line shows the average of one Monte Carlo 
run. We do not find a significant correlation between colour and overdensity, how- 
ever, there is a possible trend for a higher fraction of blue galaxies 100%) at the 
highest overdensities (log (1 + ^3) > 0.8) at all redshifts, where as the blue fraction 
at intermediate and low densities is around 80-90%. Note that there is still a 10% 
probability that this trend is caused by chance. 

Figure [3] shows the correlation between colour and local density for galaxies in 
the low and high stellar mass quartile. A clear colour offset between low- and high- 
mass galaxies of A (U — B) ~ 0.2 — 0.3 mag is present at all densities. Interestingly, 
there is a trend that mean colours of galaxies in the high-mass quartile are bluer at 



Fig. 3 Colour density relation 
in low (blue) and high (red) 
stellar mass quartiles over 
the whole redshift range 
1.5 < z < 3. The mean and 
RMS of all Monte Carlo runs 
in bins of local density are 
overplotted as red triangles 
(high quartile) and blue circles 
(low quartile). The average 
data points from the sample of 
1101 at intermediate redshift 
(0.4 < z < 1) are plotted as big 
red squares (high quartile) and 
blue pentagons (low quartile) 
respectively. 
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higher overdensities (log (1 + §3) > 0.5) than at low and average local densities. For 
galaxies in the low-mass quartile we do not see a strong correlation between colour 
and local density. 



4 Conclusions 

In this study we investigate the influence of stellar mass and local density on galaxy 
rest-frame colour at a redshift of 1.5 < z < 3. We find the following results: 

• Galaxy colour depends strongly on galaxy stellar mass at all redshifts up to z ~ 3. 
After accounting for passive evolution in colour, the colour-stellar mass relation 
does not evolve with z below that redshift. 

• Massive red galaxies exist up to z ~ 3, at the expected location of the red se- 
quence in the colour-magnitude diagram. 

• We do not find a correlation between colour and local density, however, galaxies 
in extremely high over-densities (> 5x over-dense) are bluer than galaxies in 
average and most under-dense environments. 

• More massive galaxies (logM* > 10.5M ( : ) at high relative over-densities tend 
to be on average bluer than massive galaxies at average and low local densities. 
This is due to a lack of red galaxies at high over-densities. 
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